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We introduce schemes for linear-optical quantum state generation. A quantum state generator is a device 
that prepares a desired quantum state using product inputs from photon sources, linear-optical networks, and 
post-selection using photon counters. We show that this device can be concisely described in terms of polyno- 
mial equations and unitary constraints. We illustrate the power of this language by applying the Grobner-basis 
technique along with the notion of vacuum extensions to solve the problem of how to construct a quantum 
state generator analytically for any desired state, and use methods of convex optimization to identify bounds 
to success probabilities. In particular, we disprove a conjecture concerning the preparation of the maximally 
path-entangled NOON-state by providing a counterexample using these methods, and we derive a new upper 
bound on the resources required for NOON-state generation. 
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I. INTRODUCTION 

There are many quantum states of light, which are in great 
demand in quantum technology Ufl. Due to their high ro- 
bustness to decoherence, and relatively simple manipulation 
techniques, photons are often exploited as primary carriers 
of quantum information. Recently, a number of schemes 
have been suggested enabling quantum information process- 
ing with photons using only beam splitters, phase shifters, and 
photodetectors |2|,|3j]. In a different approach to quantum com- 
puting, cluster states of photons can be used to perform com- 
putation J3l ■ Moreover, exotic states of many photons, such 
as maximally path-entangled NOON states, have found their 
place in quantum metrology |5|], lithography [6], and sensing 

A natural question arises: how can these complicated states 
of light be prepared? One approach to the problem is to make 
use of an optical nonlinearity. However, due to the relatively 
small average number of photons involved, the overall non- 
linear effect is extremely weak and typically of little practi- 
cal use (H |H QlJ • An alternative way to enable an effective 
photon-photon interaction is to use ancilla modes and projec- 
tive measurements JHHl]. In this way a quantum state gen- 
erator can be realized utilizing only linear optical elements 
(beam splitters and phase shifters) and photon counters, at the 
expense of the process becoming probabilistic ifTHl . Hybrid 
schemes, combining weak nonlinearities and measurements, 
have been also proposed recently IU2I1 . 

Despite many theoretical efforts, the problem of quantum 
state preparation of light with the help of projective measure- 
ments has not been formalized and solved in a general set- 
ting. Progress with respect to the problem of constructing 
an optimal all-optical two-qubit gate has been achieved lfl3ll . 
Moreover, a slightly different problem of attributing an effec- 
tive physical nonlinearity to a given combination of linear- 
optical transformations and projective measurements has been 
demonstrated to have a definite answer [ 14]. 

In this paper we concentrate on formalizing and solving 



the problem of quantum state generation with linear optics 
and projective measurements. We first formulate the physical 
problem of state generation in terms of polynomial equations 
in Section II. We then argue that the equations obtained can be 
solved analytically, by applying standard tools of algebraic ge- 
ometry, when a unitarity constraint is relaxed. Moreover, we 
show how unitarity can be reintroduced into the solution by 
addition of auxiliary vacuum modes. In Section III we illus- 
trate the formalism by considering an example of a "NOON" 
state generator. We demonstrate how an optimal "NOON" 
state generator can be constructed using convex optimization 
tools. A simple example for 5-photon "NOON" state genera- 
tion disproves the "No-Go" conjecture [3]. We draw conclu- 
sions and summarize in Section IV. 



II. LINEAR-OPTICAL QUANTUM STATE GENERATOR 

Any linear-optical quantum state generator (LOQSG) can 
be thought of as consisting of two main blocks (see Fig. 1). 
The first block is a TV -port linear-optical device, described by 
a unitary matrix U € U(N) combining N input into N out- 
put modes. The second block represents a projective mea- 
surement of some of the modes of the first block, in which 
a certain pattern of photons measured in some M < N of 
the modes is considered a "successful measurement", leading 
to a preparation of the desired state in the remaining modes 
(compare also Refs. HI [H H [H). This measurement is 
probabilistic — but heralded — or "event-ready". Clearly, the 
output is determined by an interplay between the numbers of 
input photons, the entries of the matrix U, and the numbers of 
photons detected. 

There are two types of problems that can be formulated 
around the concept of LOQSG. The first problem is the fol- 
lowing: Given the matrix U and a known input state, which 
output states can be generated for different projective mea- 
surements? This is what could be referred to as the "forward 
problem". This question is equivalent to the problem of find- 
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FIG. 1: The linear optical quantum state generation scheme. 



ing the effective nonlinearity generated by a given projective 
measurement and was addressed in Refs. Ill4ll7ll . The second 
problem is that of state preparation: Given an input state, a 
projective measurement, and a target state, is it possible to de- 
termine a unitary matrix U, of appropriate dimension, involv- 
ing potentially further auxiliary modes, describing the unitary 
block of the LOQSG? Given a certain input state, this impor- 
tant problem asks whether a device for the preparation of a 
certain quantum state can be identified, and if so, what ele- 
ments it contains. Once the unitary is found, there is a simple, 
well-known prescription for converting it to an optical imple- 
mentation with beam-splitters and phase-shifters, etc. flal . 
Finding such a unitary we call the "inverse problem". In this 
section, we provide a mathematical description of LOQSG 
and illustrate how methods of algebraic geometry can be used 
to solve this inverse problem. 

To be more specific, we start from a given input \ifj m ) 
in N optical modes obtained from photon sources. Then, 
we investigate all state vectors that can be reached from 
this using arbitrary networks of linear optical elements lfl8ll . 
We hence investigate the orbit SI = {\ty) : |^) = 
U(U)\ip m ) for some U 6 U(N)} of an input state vector 
\ip m ). The unitary U(U) acting in the Hilbert space of quan- 
tum states is the standard Fock-Bargmann representation of 
the unitary transformation U acting on modes. If we now de- 
note with | -0 tar ) the desired output in N-M of the modes upon 
a particular projective measurement P = \ip D ) (ip D \ on M of 
the modes, then 

9 = {|tt> : I*) = a|Ol^)+E^^r y )^f- L ).« > 0} 

i 

is the set of all state vectors which can be converted into l^' 31 ) 
by the measurement of |^i D ) in auxiliary modes. The \a\ 2 
represents a probability of success, and each \tpf ) denotes 
a state vector orthogonal to |V' D )- The solution to the state 
preparation problem is then the intersection of SI and 0. 

Since there exists a one-to-one correspondence between 
photonic states and (possibly infinite) polynomials in the cre- 
ation operators of the modes IU9I1 . the problem of finding the 
intersection of SI and can be recast in terms of polynomial 
equalities. We will assume throughout the paper that the input 
has been prepared using individual sources, so the input states 
are product states of Fock states with respect to N modes: 
\t/j m ) = \ni, . . . , njv)- These input modes are associated with 
N creation operators at = (a\, . . . , ajy). The mode trans- 
formation of the LOQSG is given by a N x N unitary ma- 
trix U. In other words, the creation operators transform as 
( a t) T = £/(a out ) T , where a out = {a\ ouV aj^J denotes 



the creation operators of the output modes. The output state 
of the LOQSG before a projective measurement hence reads 

ivn = n<o«,---,<out)io) 

n , n v m 

■■= IIs Ew i°>- « 

i=i VUl - V=i 7 

Here, F(a\ out , . . . , ajy out ) is a homogeneous polynomial of 

degree n — Ylk=i nk m ^ e crea tion operators. Coefficients 
of the monomials in F are in turn homogeneous polynomials 
of degree n belonging to the polynomial ring C[Uij]. 

The output of the LOQSG after a successful projective mea- 
surement of (mi, ■ ■ ■ , ttim) photons in modes N — M + 
1, . . . , N is then given by 

|$) = ( mi ,...,m M \r ut ). 

The state vector |$) is the result of the action of a polynomial 

^( a i oun • • • ) a ]v-A/ out) on tne vacuum |0). The polynomial 
G is constructed from F by selecting all terms which contain 

ajv-M+MuP • • • j a JV,out exact ly to me powers mi, ... , m M , 
respectively, and replacing these creation operators by 1 after- 
wards. Thus, the polynomial G acting on vacuum describes 
the final state in the first N-M unmeasured modes, which 
form the output of the device. 

For further convenience, we have summarized all symbols 
used in our paper relevant to our description of the LOQSG 
and their corresponding meaning in the table contained in Ap- 
pendix A. 

By construction G is a homogeneous polynomial of degree 
n— m, m = Y^kLi TO fc' with coefficients in C[J7,j]. Similarly, 
the desired output of the LOQSG, \ip tia ), is a state vector in 
the modes 1, . . . , N-M, and can be written as a polynomial 

Q( a \,ouv ■■■■> a ]v-M.out) actm g on vacuum. The problem of 
finding a LOQSG is then equivalent to finding a unitary matrix 
U such that G = aQ, for some nonzero a. 

We notice that for the polynomial equality G = aQ to hold 
true, the coefficients of like monomials in G and aQ must 
be equal. This leads to a system of polynomial equations 
in the variables Uij, the entries of our transformation matrix 
U. This system of polynomial equations, along with the con- 
straint that U be unitary, is the mathematical formulation of 
the state preparation problem. In general, there is no efficient 
solution technique for solving such a system. It is because the 
unitarity constraint involves the operation of complex conju- 
gation and hence cannot be written in an algebraic form, i.e., 
in the same way as the polynomial equations. So, we instead 
propose a method for finding the mode transformation matrix 
for a larger LOQSG described by input (n%, . . . , njv, 0, . . . , 0) 
and measurement (mi, . . . , tum, 0, . . . , 0). The method in- 
volves two steps: (1) For the original system, we relax the 
constraint that the matrix U be unitary, and solve the polyno- 
mial system G = aQ for non-unitary matrix A with entries 
Aij\ (2) For A as found in step (1), we find a larger matrix 
U' which is unitary and contains A as a submatrix; this U' de- 
scribes the LOQSG, and will generate the desired state |^ tar ) 
in the unmeasured modes. 
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To realize the step (1) and to solve the polynomial system 
G = aQ, we use the Grobner basis technique ll2lll . We first 
find a Grobner basis, {gi, . . . ,gi} C C[A;.j], for the ideal 
generated by the coefficients of G — aQ. This Grobner basis 
is the minimal set of polynomials that do not leave a remain- 
der when the original polynomials are divided by them. One 
of the essential features of this approach is contained in the 
elimination theorem, which states that the Grobner basis con- 
sists of polynomials of only the first variable in the ordering, 

gi=gi(A ltl ) cC[Ai,i], 

the first two variables, g 2 = 92(^1,1)^1,2) C CfA^i, Ai 2], 
and so forth. This property allows us to simply find the 
solution to the whole system by solving for each variable 
subsequently — a procedure that is often regarded as exten- 
sion of a solution. That means finding the solutions A\ 1 of 
9l (A M ) = 0, A\ 2 of 32(^1,1,^1,2) =0, and so on. Thus, 
the whole problem is reduced to finding roots of monovari- 
ate polynomials. In this way, the Grobner basis technique is 
the generalization to polynomial systems of the technique of 
Gaussian elimination for linear systems, and is simple to im- 
plement in algebraic software. 

For the case of product-state inputs, we can simplify the 
polynomial system of equations by eliminating some arbitrary 
factors before using the Grobner basis technique. We remark 
that for any matrix A, which is a solution to the polynomial 
system, a rescaling of the rows of A — corresponding to the 
input modes of the product state — by arbitrary nonzero fac- 
tors Xi and of the columns of A — corresponding to the mea- 
sured modes — by arbitrary nonzero factors yj will gener- 
ate another solution to the system. Hence, for any solution 
A, we define the equivalence class of A by [A] = {XAY : 
X = diag(x!, . . . , x N ), Y = diag(y x , . . . , y N ),Vi = ■■■ = 
Un-m = !}■ In any such equivalence class, there is ex- 
actly one member A for which A11 = ... = Ajy^ = 
Ai,jv-a/+i = . . . = Ai,jv = 1. Hence, by setting the 
aforementioned variables to 1 from the very beginning, one 
can solve a much simpler polynomial system in the remaining 
variables using the Grobner basis technique; each solution of 
that system will correspond to an equivalence class of solu- 
tions for the entire system. We call such a solution an equiva- 
lence class representative. 

For a given input and projective measurement, one will find, 
in general, many equivalence classes of solutions. In most 
cases, an equivalence class representative A will not be uni- 
tary and hence will not conserve the canonical commutation 
relations, which implies that A and hence an underlying LO- 
QSG is not physically realizable. However, for a non-unitary 
A, one can think of a larger optical system described by a uni- 
tary matrix U', and hence experimentally implementable, in 
N' = N + d modes which in turn contains a member of [A] 
as a submatrix. Physically, this corresponds to an addition of 
auxiliary modes to the LOQSG. Then, if we input and mea- 
sure vacuum in these additional modes, the entries of U cor- 
responding to the added modes do not affect the final state of 
the original modes. Given an input on N modes and a desired 
output on N — M modes, the (i) total number of modes N', 
the (ii) unitary U' S U (N 1 ), and the (iii) measurement pattern 



(mi, . . . , tum-, 0, . . . , 0) reflecting success are then the solu- 
tion to the problem of finding a LOQSG. The question arises: 
for a given non-unitary matrix A, when is it possible to find 
a matrix U' which is unitary and contains some member of 
[A] as a submatrix? Moreover, how can we find the extension 
which optimizes the success probability of the LOQSG? 

It is not difficult to see that any matrix A can be extended 
if A = 1 1 1 A| 1 1 < 1, so if the largest singular value of A is not 
larger than unity (compare also Ref. IU3IP . Hence for any A, 
we can always extend A/X 6 [A] to a unitary. Let the singular 
value decomposition of A be denoted as A = VDW, where 
V, W € U(N) are unitary and D = diag(di, . . . , (In). Then, 
the minimum dimension of the extended unitary U containing 
A/X is N + d, where d = rank(L>/A - 1). Thus, d < N 
additional vacuum modes are sufficient to extend a member of 
[A] to a unitary. The minimum number of additional modes 
needed to extend at least one member of [A] remains unsolved. 
However, we can solve the problem of optimizing the success 
probability of the extended LOQSG using, for example, ideas 
from optimization theory 12311 . 

For any solution of the polynomial equations the success 
probability of correct state preparation is by definition given 
by \a\ 2 . The allowed rescaling of rows and columns rescale 
the amplitudes in the respective input and output modes such 
that within an equivalence class 
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where p s is the success probability and Xk and yi are the arbi- 
trary row and column multipliers as presented in the definition 
of our equivalence classes. Note, that w.l.o.g. it is sufficient 
to chose Xk and yi real and positive, thus every quantity in 
(O is real. The problem is then to maximize p s subject to the 
constraint X AY 2 A^ X < 1. For simplicity we will focus on 
the important case when m — rij for alH, j = 1, . . . , N and 
rrii = rrij for all i, j = 1, . . . , M. Then, for given Y, the con- 
straint can be written as a so-called semi-definite constraint. 
Semi-definite optimization problems can be efficiently solved 
and solvers are readily available 12411 . As AY 2 A' is invertible, 
the constraint is equivalent to 



1 X 

X (AY 2 A^)- 1 



> 0. 



The objective function, p s = c(J\ k Xk) f° r a constant c 



2 (IW 



is a monomial, and hence clearly not linear. 



However, this monomial can be relaxed to again a hierarchy 
of semi-definite constraints, without altering the optimal ob- 
jective value: Let s be the smallest integer such that 2 s > N. 
Let us assume that 2 s = N; if N is smaller, we can always 
pad with variables that we enforce to be unity by means of 



linear constraints. Then, let x k = x k for k = 1, 



.N, 



and x 



U) Jj) 

k-l x k 



= (4/2 ) 2 forfc e {2,4, ...,2 S+1 -J'} and 
j = 1, . . . , s. Hence, we have introduced a number of new 
variables, according to a hierarchy. Each of the quadratic 
equality constraints of this form can actually be enforced, as 
is easy to see, see footnote 12511 . In fact, given Y, we have 
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written the above problem as a semi-definite problem. In turn, 
relaxations give efficient upper bounds to the success proba- 
bility for the LOQSG when simultaneously varying X and Y. 
The solution to this problem gives rise to the LOQSG operat- 
ing with the highest probability of success within the equiva- 
lence class of A. Then, for a given input and measurement, the 
overall optimal LOQSG can be found by simply optimizing 
within each equivalence class by the method described above, 
and choosing the best equivalence class. 



III. EXAMPLE: NOON-STATE LOQSG 

To illustrate how to solve the problem of the intersection of 
and 9 using the language of polynomials, let us consider 
the state generation problem for target states of the form 



and G. They read, 



\^} = (\n,0) + \0,n))/V2, 



the so-called path-entangled NOON state [20J]. The NOON 
state is important in a number of applications such as quan- 
tum metrology and lithography [6]. It cannot, however, be 
generated using only linear optics from product sources, and, 
therefore, the NOON-state generation problem is both impor- 
tant and nontrivial. For NOON-state generation, we require 



F(Kout) n + (4,out)")- 



TV 

Coefficient^, n(«Lt) m3 ] = - 

j=3 V 

(3) 

Eq. (f3]i cannot be satisfied for input states for which rii > 
rn + 1 for any i E {1, . . . , N}. To see this, note that the RHS 
of Eq. (O has a unique factorization into distinct linear terms. 
Ifrii > m + 1 for any mode i the LHS would necessarily have 
multiple identical linear factors of in a\ out and a\ out . Hence, 
for the equality to hold, we must have rii < m + 1. Thus, 
using polynomial properties, we have shown that the intersec- 
tion of fi and 6 is empty for any \ip m ) which contains more 
than TO+l photons in any of its modes. In other words, we 
have demonstrated that the largest NOON state that might be 
generated from N modes with a total of m measured photons 
is N(m + 1) — to. 

The argument above is only useful for determining whether 
an intersection is empty. However, the ultimate goal of this 
section is to present a constructive method to find all points 
in the intersection when it is non-empty. To illustrate our 
method, we demonstrate how to build a NOON-state LOQSG. 
Let us consider the following LOQSG input \ip in ) = |2, 2, 2). 
In this case the mode transformation matrix A is a 3 x 3 ma- 
trix with 9 complex entries Aij. We are particularly inter- 
ested in a projective measurement of one photon in mode 3, 
i.e., \ip D ) — \mi) = The total number of photons in the 
remaining (unmeasured) modes is thus 5 and, therefore, the 
only NOON state that can be created is the one corresponding 
to the state vector \ip Vdr ) = (|5,0) + |0, 5))/V2. Proceeding 
as discussed in Section II we calculate the polynomials F, Q, 



G = Coefficient[F(a| out ,4 01Jt ,a 

" -J= ((4 >out ) 5 + (4 )OUt )5 



4, out)' 4, out]' (4) 



V2~5! 

The equality G = aQ leads to a system of six polynomial 
equations for the coefficients of (o\ out ) k (a\ out ) 5 ~ k , k = 
0, . . . , 5, with respect to ten complex variables: the nine el- 
ements of A and a. 

We now simplify the system to solve only for equivalence 
classes by setting Ai i = A-2,i = A31 = A13 = 1. Once 
the Grobner basis technique is used to solve the simplified sys- 
tem, one finds that there are finitely many equivalence classes. 
Here, we will illustratively show the work for one of them (in 
fact, for an optimal one). The following matrix is the optimal 
equivalence class representative: 



A = 



1 



e 5 



1 

3+vy 
2 e 

ih/5, 
2 c 



(5) 



Although the mode transformation in Eq. (0 generates the 
final desired NOON state vector up to normalization 

a(|5,0) + |0,5))/V2 

with a = ^30(3 + v5), it is easy to see that the matrix A is 
non-unitary and must be extended. 

We now seek to optimize the success probability, p s = 
a 2 (xiX2X3) 4 (y3) 2 of the extended unitary U over all 
Xi,X2,%3, and 7/3 such that XAY is extendable. Let 
us view xi,X2,xz in spherical coordinates: p s — 
a 2 r 12 (sin 2 9 cos 9 sin <fi cos </>) 4 ?/|. For a given direction 9, </>, 
the optimal p s is reached when r is chosen maximally; thus, 
we choose r so that the largest singular value of XAY is 1 . 
For a given y 3 , then, the optimal success probability of the 
extension over all X can be found by varying 9 and <f> over 
an entire sphere; call this p s .opt- Then, by evaluating p Si0p t 
over a reasonably large range of y 3 , the global optimum of the 
extension can be found. 

Using the optimization procedure described above, we find 
that the global maximal success probability is reached at a 
point [x\, X2, X3, 2/3) for which the matrix can be extended 
minimally, i.e., by only 1 additional dimension. It remains 
unclear whether this result will hold true in general. One can 
check that the unitary matrix 



0.5722 0.5722 

0.5257 0.5257c- 1 

0.5257 0.5257^ 

0.3461e" 0.34616^ 



0.1894 
0.45566"^ 
0.4556e¥ 

0.7409 



0.5561 
0.4895e^ 
0.4895e-- 

0.4599 



contains a member of [A] as a submatrix and generates a uni- 
tary transformation on four modes which results in a five- 
photon NOON state starting with input 1 2,2,2,0) and pro- 
jective measurement of one photon in mode 3 and vacuum in 
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mode 4 with optimal success probability sa 0.05639. By using 
the well-known algorithm presented in Ref. [18], the above 
unitary matrix can easily be transformed into a linear-optical 
network consisting of phase-shifters and beam-splitters and, 
as such, can be constructed in the lab. 

Remarkably, the example we have just considered also 
serves as a counterexample to the "No-Go" conjecture @]. 
The conjecture states that the largest NOON state that can be 
produced from product inputs with only TV modes is that of TV 
photons. However, we have seen that there is a situation when 
the NOON state of live photons can be generated in only four 
modes. 



IV. CONCLUSIONS AND OUTLOOK 

In this work, we have shown how the problem of identifying 
a linear optical state preparation device can be formulated and 
solved using the language of polynomials. In this way, we do 
not have to include unitarity conditions as constraints on our 
system. Instead, we solve the polynomial equations using the 
methods of algebraic geometry and later restore unitarity with 
a sacrifice of at most doubling the number of modes. We have 
introduced a general framework that allows for the system- 
atic construction of linear optical devices preparing entangled 
quantum states of light. 

It is worth noting that in addition to solving the problem 
of state generation by linear-optical means, the technique we 
propose can also be used to construct optimal linear-optical 
quantum gates. With care taken to formulate appropriate 
equivalence classes, the problem of constructing a (probabilis- 
tic) linear-optical quantum gate involves the same system of 
equations, and hence can be solved using Grobner bases and 
vacuum extensions. Applied to gate construction, our solu- 
tion scheme can be seen as the generalization of the proce- 
dure used to construct an optimal linear-optical NS gate in 
Ref. 01. 



V. APPENDIX A 



Symbol 



N 
M 
U 
TH 
n 

rrii 
m 

l^ m ) 
|i/> out ) 

m 

\m 

F 
G 

Q 
a 
A 
d 

N' 
U' 



Meaning 



total number of modes 

number of ancillary modes 

TV-dim. unitary representing the linear optics 

number of photons input in mode i 



total number of input photons 



number of photons measured in mode TV — M 



i = N-M+l 



- detected photons number 



TV-dim. input product state 

TV-dim. state after linear optics 

M -dim. specified measurement state 

TV — M-dim. state remaining after measurement 

TV — M-dim. desired output state 

polynomial representing |?/; out ) 

polynomial representing |$) 

polynomial representing I?// 31 ) 

probability amplitude for successful measurement 

TV-dim. matrix solution to G = aQ 

number of vacuum modes added to system 

N + d 

TV'-dim. unitary containing A as a submatrix 
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